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Abstract 

We demonstrate the theoretical possibihty of obtaining a pure spin current in a ID ring with 
spin-orbit interaction by irradiation with a non-adiabatic, two-component terahertz laser pulse, 
whose spatial asymmetry is reflected by an internal dephasing angle (j). The stationary solutions of 
the equation of motion for the density operator are obtained for a spin-orbit coupling linear in the 
electron momentum (Rashba) and used to calculate the time-dependent charge and spin currents. 
We flnd that there are critical values of (j) at which the charge current disappears, while the spin 
current reaches a maximum or a minimum value. 
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I. INTRODUCTION 



Obtaining and controlling spin currents in solid structures has been one of the most im- 
portant goals of spintronics research. In the recent past, the focus of this endeavor has been 
on the creative use of the spin-orbit interaction (SOI) that appears in systems with broken 
inversion symmetry, be that via confinement (Rashba)Q or in the bulk (Dresselhaus) j2|. The 
well known properties of quasi-one dimensional rings to support persistent charge currents 
[sl-lsj made them particularly appealing for spin-dependent exploration. The underlying 
physical phenomenon leading to the creation of persistent currents is the imbalance in the 
left/right charge carrying states realized, initially, in the presence of a magnetic flux threaded 
through the center of the ring. In a similar experimental setup, the existence of persistent 
spin-currents was obtained in the presence of SOI 0-0], when electron-electron interaction 
are neglected. Both analytical joj and numerical results jsf, point out that in the presence 
of a magnetic flux, the charge and spin currents in the rings are simultaneously present. 

In this paper, we revisit this problem from the perspective of creating persistent charge 
and spin current by non-adiabatic methods. As before, the original exploration of such 
ideas was focused on the generation of charge currents various mechanisms able to realize 
of an imbalance between leftright momentum-carrying states [o - 11 1. In particular, we are 



concerned with the application of a ultrashort, terahertz frequency laser pulse, endowed 
with a spatial asymmetry expressed through an internal dephasing angle. The efficiency of 



this method for creating persistent charge currents has been explored in Refs. [Ill, Il2|. Its 
extension to a ring endowed with SOI is discussed below. As we demonstrate, the interplay 
between the spin orbit coupling that rotates the electron-spin around the ring and the spatial 
asymmetry of the initial excitation generates favorable conditions where the charge current 
disappears, while the spin current reaches a maximum level. For both a dipolar and a 
quadrupolar perturbation, we can find the critical values of the dephasing angle for which 
this situation occurs. 



II. THE RING MODEL 

We consider a one-dimensional (ID) quantum ring of radius tq containing few electrons, 
endowed with a Rashba interaction, linear in the electron momentum. In a discrete (tight- 
binding) representation the ring is reduced to N sites (points) distributed on a circle, whose 
angular coordinate is given by On = 2mT /N with n = 1, ...,N the site index. The Hamiltonian 
describing the noninteracting electrons is written in terms of the creation and annihilation 
operators c^^ and associated with the single-particle states \na), where a = ±1 is the 
spin index. This Hamiltonian has been extensively discussed in literature [6|-@, 13|, so here 
we will write it directly: 

n,a n,cr J 

- iVa J2 K(^n,n+l)],,, 4,C„+i^, -f-h.C. (l) 
n,a,a' 

The two energy scales of the problem are set by the hopping matrix element V = 
/2m*a^, where m* is the effective electron mass and a = 27rro/N is the discretization 
constant, and by the Rashba coupling, of strength a, which generates Va = a/2a. The spin 
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operator ar{0) introduced in Eq. ([T]) represents the local orientation of the electron spin 
along the radius of the ring and is given by a linear combination of the Pauli operators 
written for the azimuthal coordinate 9n,n+i = {On + dn+i)/'2, 

(^r{On,n+l) = CTx COS dn,n+l + CTy siu dn,n+l ■ (2) 

The energy spectrum of the Hamiltonian, calculated for an even number of sites, N = 20, is 




FIG. 1: (Color online) The eigenvalues Ei± of the ID quantum ring with N=20 points vs. quantum 
number /. (a) Va = 0.0, (b) Va = 1.0. The energy unit is V = h'^/2m*a'^. 



shown in Fig. [T]for Va = and for Va = 1.0 units of V. The split realized by SOI is described 
by two eigenvalues, (right arm) and Ei_ (left arm), with / = 0, ±1, ±2, ■ ■ ■ , ±{N/2 — 
I), N/2, 

El + Ei±i El — Ei±i 



Eli = ' ^ '"^ + ' ^ '"V l + tan^ 29a , (3) 

where Ei = 2V — 2^ cos(27r//A^) is the degenerate eigenvalue in the absence of SOI and 9a 
is given by 

The corresponding eigenvectors of the Hamiltonian ([T]) are: 

i*H>^^i:e-(_;r,t,Ji">. 



IT \ ^ iif) f e *^"sin6'„\ I , 



The velocity operator is given by 
vg = r9 = r^[H,9] 



--ri'^Y. cLc„+la + 17^ K(^n,n+l)]<,,, 4<xC„+l.' + h.C.} , (6) 
It n,a ^ ' 
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FIG. 2: Spectrum of the ID quantum ring with N=20 points vs. Rashba parameter V^. The 
eigenvalues and Eq± are dotted to illustrate the lifting of degeneracy when Va / 0. 

and its eigenvectors are also |^I^i±), with the associated eigenvalues 



vi± = \ yl + tan^ 26^ , (7) 

where vi = 2 (Va/h) sin(27r//A^) is the velocity in the absence of SOI. 
It is customary to introduce the tilt-spin operator 829^, 

S29c = COS 29aS^ - sin 26'«S'r, (8) 

which is a linear combination of Sz and Sr, the spin operators for z and radial directions, 
respectively. The spinors \'^i±) are also eigenvectors of 829^ associated with eigenvalues 

±h/2. 

In the absence of SOI, all quantum states are four-fold degenerate, except those at the 
edges of the spectrum which are only twice degenerate. This can be seen in Fig. [U^a) where 
for Va = the two branches Ei+ and Ei_ coincide. In this case the spin operator 829^ 
becomes 8z and \'^i±) become the |t) and \l) eigenstates of 8z. For 7^ 0, the energy 
spectrum becomes broader and all states are twice degenerate since Ei± = E^i^. States 
at the crossing points in the spectrum shown in Fig. [2] preserve the four-fold degeneracy. 
Since these degenerate states carry opposite currents vi±. = — W-/^, the system can support 
a nonzero spin current. 



III. THE TIME EVOLUTION 

At t = the quantum ring described above is exposed to a short terahertz two-component 
pulse, 

Hn{t) = Ae~^^ [sin(wit) cos6' + sm{uj2t) cosn{6 + 0)] , (9) 

of duration ~ F^^ and amplitude A [Hj. n = 1,2 describes the multipole order of the 
second component, while the dephasing angle between the two components makes the ex- 
ternal perturbation asymmetric. The terahertz scale of the excitation frequencies 071^2 is at 
least an order of magnitude larger than the spin relaxation rates in InAs semiconductor het- 
erostructures which reaches values from anywhere between tens and hundreds of picoseconds 



14, 15 . 
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The time evolution of the system's observables is determined by using the density operator 
p{t), which at t > satisfies a quantum Liouville equation 

thp{t) = [H + Hn{t), p{t)] . (10) 

The initial condition is such that p{t = 0) represents the ground-state density operator 
which is constructed in terms of the eigenvectors | "^i^j) shown in Eq. of the initial, 
time- independent Hamiltonian, Eq. ([1]): 

p{t = 0)=Y.Pl<r\'^la){^la\, (11) 

1(7 

where pia^ are the populations of the degenerate single-particle states {J2iaPicT = !)• 

For any t > 0, Eq. fITU]) is solved numerically and p{t) is obtained by using the Crank- 
Nicolson finite difference method [lH with small time steps 6t <^ F"^/. The expectation 
value of any observable O is then calculated as (O) = Tr {p{t)0). 



IV. RESULTS 

We define the charge current as I'^ = evg and the spin current along a direction as 
K = ^ {a^ve + vgau). Since the spin operator S2ea& commutes with the unperturbed 
Hamiltonian ([1]), we calculate the expectation value of which corresponds to the direction 
of the spin e2g^ = cos2^Q,e^ — sin26'Q,er. To simplify the notation, we denote /*(t) = (lig^), 
vg(t) = {vg) and I'^{t) = evg{t). Since vg and I^g^ commute with H, the charge and spin 
currents become constant after the external perturbation vanishes. 

To illustrate our results, we consider an InAs (electron effective mass is m* = 0.023me) 
quantum ring of radius tq = 14 nm. We choose our Rashba parameter Va = 0.05, which 
corresponds to a SOI strength a = 37.56 meVnm, within the range of experimentally de- 



termined values [16|. The number of sites used in the discretization is = 20, leading to a 
length unit a = 4.4 nm and an energy unit V = 85.6 meV. For a system with Ue = Q electrons, 
in the many-particle, non-interacting ground-state, the occupied states are '^o±, '^i± and 
^-i±5 with weights po± = Pi± = P-i± = l/''^e (and all the other pi^ = 0). The many-particle 
ground-state is therefore nondegenerate. In this configuration, i.e. at t = 0, the average ve- 
locity is vg(t = 0) = and the spin current P(t = 0) = 2/1(^0+ + + f-i+) = ~0.089Va/h. 

The first external pulse we consider is the superposition of two dipoles, corresponding to 
77, = 1 in Eq. Q. For the selected parameters of the ring, in the absence of the SOI, we obtain 
the Bohr frequencies huo^i = 2.89 meV and huji^2 = 8.60 meV (defined as ui^i> = \Ei — Eii\/K). 
In Figs. [3] and |4] we show the numerical results obtained for frequencies huji = 2.83 meV, 
hu2 = 8.11 meV, with the attenuation factor F = 4ui, and amplitude A = 67.68 meV. The 
duration of the pulse is t/ ~ 0.5 ps. Since the pulse produces many-particle excited states, 
in our calculations uji and uj2 are chosen to be slightly different from the Bohr frequencies, 
in an effort to create a more realistic algorithm that reproduces closely what happens in an 
experimental situation. Moreover, for ui and U2 exactly equal to the Bohr frequencies, the 
outcome is preserved. 

The time evolution of velocity vg{t), proportional with the charge current I'^{t), and that 
of the spin current P{t) are illustrated in Fig. [3] for dephasing angles = 0, 7r/2 and vr. In 



Fig.[3](a), we recover the result of Ref.llll where it was demonstrated that a charge current can 
be non-adiabatically generated through the application of a spatially asymmetric terahertz 
excitation. Thus, for (p = 0, vg = 0, while for = n/2, vg ~ 0.37Va/h. 
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FIG. 3: (a) The time evolution of the azimuthal velocity v(t) = (vg) = I^{t)/e, (b) the spin current 
P{t) = (120^), for the selected dephasing angles 4> = 0, 7r/2, vr. The inset shows the two components 
of the radiation pulse described by Eq. ([9]) with n = 1, for 6 = (p = 0. The lower peak corresponds 
to the first term and the higher peak to the second term. 



In Figure |3]^b) we present the time evolution of the spin current, /^(t), corresponding to 
the spin projected on the proper axis 626*^. As previously stated, in the presence of SOI, the 
initial state of the system has a nonzero spin current, /''(O) ^ 0. On account of the external 
pulse, the electrons in the ring are non-adiabatically excited, and the spin current evolves 
towards a new steady-state value, = I^{t ^ tj). After the external pulse vanishes, the 
spin current Is{t) is constant in time, but its amplitude varies with the dephasing angle. 

We call vg and l2e the constant values of the velocity and spin current after the perturba- 
tion vanished (i.e. vg(t ^ t/) and P(t S> t/) respectively). They are plotted as a function of 
the dephasing angle </> G [0, 27r] in Fig. HJ but we note that they reflect the periodicity in (p of 
the applied pulse. For (p^i = O.Gtt and (j)v2 = 1.47r the charge current is maximum and mini- 
mum, respectively. For equal to and vr no charge current flows through the ring [vg = 0) 
whereas a spin current is present, I20 7^ 0, being maximum for (psi = and minimum for 
(f)s2 = TT. These are the critical dephasing angles for which only pure spin currents exist. In 
this case, the numerical results show that the states with opposite velocity, "^i^ and \E'(_;)_, 
are equally excited by the external pulse netting zero total velocity and no charge current 
in the steady state. This is not the case for intermediate angles where states with opposite 
velocity are asymmetrically excited and consequently, generating nonzero values for both 
charge and spin currents. 

In Fig. [5] we display vg and l2e obtained by exciting the system with a pulse given by 
Eq. (Q written for n = 2, i.e. by a combination of a dipole and phase shifted quadrupole. 
The frequencies ui, U2, the attenuation factor F, and the amplitude A remain the same 
as before. Due to the quadrupolar component of the pulse, the period of vg and l2g is 
halved to A0 = tt. The maximum and minimum values of vg are reached for = 0.287r 
and (f)y2 = 0.727r (nearly half of the previous values). For the critical angles 0^1 = and 
(ps2 = 7r/2 no charge current is induced in the ring {vg = 0), but the induced spin current 
reaches extreme values, minima or maxima respectively. Again, in these situations only a 
pure spin current is induced in the ring. 
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FIG. 4: The stationary values of the velocity ve and spin current /|g after the external pulse 
vanishes (at t 3> t/) vs. the dephasing angle (p. The external pulse is Hi{t), Eq. ([9]) with n = 1. 
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FIG. 5: The stationary values of the velocity vq and spin current /|g after the external pulse 
vanishes (at t ^ tf) vs. the dephasing angle (j). The external pulse is H2{t), Eq. Q with n = 2. 



V. CONCLUSIONS 

In conclusion, we studied the non-adiabatic excitation of spin and charge currents in a 
ID quantum ring in the presence of Rashba SOI. The ring was subjected to an external 
pulse that is spatially asymmetric, having two components with a relative dephasing angle. 
We investigated two models, a dipole plus a rotated dipole and a dipole plus a rotated 
quadrupole. By numerical calculation, we showed that for certain values of </> called 0„i and 
0^2 the induced charge current reaches extreme values with P{(f)yi) = —P{(f)v2)- Due to the 
presence of SOI a nonzero spin current is also induced in the ring, for both pulse models, with 
amplitudes depending on the parameters of the pulse. We found simple critical values of the 
dephasing angle, 0ci and 0c2, for which the induced charge current disappears, whereas the 
spin current reaches maxima or minima. The method may be used in practice to convert 
an optical signal of variable amplitude into a dissipationless (persistent) spin current for 
information transfer purposes. Similar results can be obtained with the Dresselhaus SOI 
instead of the Rashba SOI, due to the equivalence of the corresponding Hamiltonians. 
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